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ABSTRACT 

A puzzling feature of the Chandra - detected quasar jets is that their X-ray emission decreases faster along the 
jet than their radio emission, resulting /Shan outward-increasing radio-to-X-ray ratio. In some sources this behavior 
is so extreme that the radio emission peak is located clearly downstream of that of the X-rays. This is a rather 
unanticipated behavior given that the inverse Compton nature of the X-rays and the synchrotron radio emission 
are attributed to roughly the same electrons of the jet’s nonthermal electron distribution. In this Letter we show 
that this morphological behavior can result from the gradual deceleration of a relativistic flow and that the offsets 
in peak emission at different wavelengths carry the imprint of this deceleration. This notion is consistent with 
another recent finding, namely, that the jets feeding the terminal hot spots of powerful radio galaxies and quasars 
are still relativistic with Lorentz factors T ~ 2-3. The picture of the kinematics of powerful jets emerging from 
these considerations is that they remain relativistic as they gradually decelerate from kiloparsec scales to the hot 
spots, where, in a final collision with the intergalactic medium, they slow down rapidly to the subrelativistic 
qi velocities of the hot spot advance speed. 

Subject headings: galaxies: active — quasars: general — radiation mechanisms: nonthermal — X-rays: galaxies 


1 . INTRODUCTION 

The superior angular resolution and sensitivity of Chandra 
has led to the discovery of X-ray emission from a number of 
quasar jets. Schwartz et al. (2000) were the first to note that 
the X-ray emission from the knots of the jet of the superluminal 
q2 quasar PKS 0637-752 (Chaitas et al. 2000), at a projected 
distance ~100 kpc from the quasar core, is part of a spectral 
component that is separate from its synchrotron radio-optical 
emission and is too bright to be explained by synchrotron self- 
Compton emission from electrons in energy equipartition with 
the jet magnetic field. Note, however, that in the innermost knot 
of some sources — e.g., 3C 273 (Marshall et al. 2001) and PKS 
1136-165 (Sambruna et al. 2002) — a synchrotron X-ray con- 
tribution is possible. These properties are apparently common 
to other quasars jets, as indicated by the mounting observational 
evidence (Sambruna et al. 2001; Marshall et al. 2001; Jester 
et al. 2002; Sambruna et al. 2002; Siemiginowska et al. 2002; 
Jorstad, Marscher, & McHardy 2003; Siemiginowska et al. 
2003; Yuan et al. 2003; Cheung 2004; Sambruna et al. 2004). 

To account for the level of the observed X-ray emission, 
Tavecchio et al. (2000) and Celotti, Ghisellini, & Chiaberge 
(2001) proposed that this is due to external Compton (EC) 
scattering of the cosmic microwave background (CMB) pho- 
tons off relativistic electrons in the jet, provided that the jet 
flow is sufficiently relativistic (T - 10) to boost the CMB en- 
ergy density in the flow frame (by T 2 ) at the level needed to 
reproduce the observed X-ray flux. 

In all these sources the radio-to-X-ray ratio increases down- 
stream along the jet, an unexpected behavior, given that the 
cooling length of the EC X-ray-emitting electrons (7 of about 
a few hundred) is longer than that of the radio-emitting ones 
q3 (7 of about a few thousand) and comparable to or longer than 
the size of the jet, which would lead to a constant X-ray bright- 
ness as far out as the hot spots, contrary to observations. More 
surprisingly, in some jets (e.g., 3C 273 in Sambruna et al. 2001 
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and Marshall et al. 2001; PKS 1136-135 and 1354+195 in 
Sambruna et al. 2002; PKS 1 127 — 145 in Siemiginowska et al. 
2002; and 0827+243 in Jorstad et al. 2003) the X-ray and radio 
maps are anticorrelated, with the X-ray emission peaking closer 
to die core, gradually decreasing outward, and the radio emis- 
sion increasing outward to peak practically at the very end of 
the radio jet. 

To explain the reduction of the X-ray flux along the jet, 
Tavecchio, Ghisellini, & Celotti (2003) suggested that the X- 
ray emission originates in a collection of microknots under- 
going adiabatic expansion sufficient to produce the desirable 
electron cooling. However, this would also suppress the radio 
emission, leading to practically indistinguishable radio and X- 
ray morphologies, contrary to their observed spatial anticor- 
relation. An elegant suggestion by Dermer & Atoyan (2002) 
that the X-rays are synchrotron emission from electrons cooling 
in the Klein-Nishina regime naturally produces shorter sizes in 
X-rays than in radio. However, it also produces larger optical 
than X-ray sizes, contrary to observations, and it seems there- 
fore not to be applicable in this particular context. 

Our view is that although X-ray-producing electrons are 
present throughout the jet, the X-ray brightness decreases be- 
cause the CMB photon energy density in the flow frame de- 
creases along the jet as a result of a decelerating relativistic jet 
flow. The decrease in the flow Lorentz factor leads to a decrease 
in the comoving CMB photon energy density and hence to a 
decrease in the X-ray brightness along the jet. At the same 
time, the flow deceleration leads to a compression of the mag- 
netic field; this results in an enhanced radio emission with 
distance, which thus gets displaced downstream of the EC X- 
ray emission. 

Based on the radio, optical, and X-ray jet maps, we argue 
in this Letter that powerful extragalactic jets are relativistic and 
gradually decelerating. In § 2 we formulate the synchrotron 
and EC emission process from a decelerating jet flow and pre- 
sent our results. In § 3 we discuss our findings and touch upon 
some open issues. 
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Fig. 1. — EC ( upper panel ) and synchrotron (lower panel ) flux normalized 
to its maximum value for a range of observing angles along a jet that decelerates 
from r o = 6 to T = 2 with e = 1 from = 100 kpc to = 300 kpc, while 
its radius doubles. The observing angles are 0 = 0° ( short-dashed line), 
6 = 1/2 T 0 w 4?8 ( long-dashed line), 6 = 1/T 0 9°6 (solid line), and 6 = 

2/r o =s 19°5 (dot-dashed line). 

2. THE MODEL 

We parameterize the jet flow assuming that its radius r scales 
as r{z) - r 0 /(z), where /(z) is a function with f(z 0 ) = 1 and 
r 0 is the jet radius at a fiducial point z 0 along the jet. Similarly, 
we assume that the flow decelerates as T(z) = T 0 g(z), where 
g(z) is a function with g(z 0 ) ~ 1 anci T 0 is the bulk Lorentz 
factor at z 0 - The evolution of the electron energy distribution 
(EED) along the jet is determined by the combination of ad- 
iabatic changes (losses for expansion and gains for compres- 
sion), radiative losses, and particle acceleration. If the radiative 
losses of the electrons responsible for the X-ray and radio 
emission are negligible (as may be the case for the (7 of about 
a few hundred) EC X-ray-emitting electrons, or if we assume 
that some localized (e.g., shocks; Kirk et al. 2000) or spatially 
distributed (e.g., turbulence; Manolakou, Anastasiadis, & Vla- 
hos 1999) particle acceleration mechanism offsets the radiative 
losses, then the evolution of the EED can be treated as adiabatic. 
In this case, assuming that T » 1, the elementary volume dV 
scales as dV = dV 0 f 2 g , allowing for expansion or compres- 
sion, while the adiabatic electron energy change rate with z is 
obtained from 7 ' = -7 (2 f'ff+ g'lg ) /3, where the prime de- 
notes differentiation with respect to z. The solution of the above 
equation is 7 (z) = y 0 f~ z/3 g ~ m • Assuming particle conserva- 
tion, n(y, z)dV dy = n(y 0 , z 0 )dV 0 dy 0 , the comoving EED 
n(% z) can be written as 

n(7,z) = *7~ 5 r 2( ' +2) V<' +2V3 , (1) 

where 77(7, z 0 ) = ky~ 5 is the EED at z 0 - Following Georgan- 
opoulos, Kirk, & Mastichiadis (2001), the EED in the local 
rest frame is n( 7, z)A 2+s , where A(z) is the familiar Doppler 
factor A(z) = 1/ (T(z) [1 — (3(z ) cos 0]} and 6 is the observing 
angle. Using this and taking into account the cosmological 
corrections, the EC flux per dz in the 5-function approximation 
(Coppi & Blandford 1990; an electron of Lorentz factor 7 


upscatters seed photons of energy e 0 only to an energy 
4e o7 2 /3) is 

dK_ = ko r cU{ 1 +Z) 2 ,/3e\- ( - w 

dt dz 2 d% J S \4e 0 J 

(2) 

where Z and d are, respectively the redshift and proper distance 
of the source, € 0 and e are, respectively, the energy of the CMB 
seed photons and the observed photons in units of m e c 2 at 
Z = 0, and U is the CMB energy density at Z = 0. For 
0 = 0, A (z) = 2r(z), and the z-dependence of the flux scales 
proportional to 4+2 ' y3 . 

Assuming flux conservation, the magnetic field B{z) can be 
written as B(z) = B 0 /fg. The synchrotron flux per dz in the 5- 
function approximation (an electron of Lorentz factor 7 pro- 
duces synchrotron photons only at the critical synchrotron en- 
ergy By 2 IB CT ) is then 

dF s _ ka T cB cr Br iy2 
de dz 12tt^ 2 (1 + Z ) ( * +1)/2 

x y--(75— l)/6g— <5* + 7)/6^ + 3)/2^ e £ y(s- l)/2 ^) 

where B ct = m 2 c 3 leh = 4.414 x 10 13 G is the critical mag- 
netic field. Note that while the synchrotron flux decreases with 
increasing redshift [the quantity d 2 { 1 + Z) (5+1)/2 increases], the 
EC flux [oc (1 + Z) 2 ld 2 ], after an initial decrease up to Z» 

1.9, remains practically constant for higher Z, as was first noted 
by Schwartz (2002). 

For 6 = 0 the z-dependence of the flux scales proportional q 4 
to f-Vt-Wg-i *- 1 > /3 . The condition under which the synchroton 
flux increases along the jet at 6 = 0 is that deceleration dom- 
inates over the jet opening, leading to -g'lg < (7s - 1) x 
f'!2(s - 1)/. Taking the ratio of equation (3) to equation (2), 
the condition under which L S IL BC (or, equivalently, the radio- 
to-X-ray spectral index a rx ) increases with z at 6 = 0 is 
—g'lg and for a given jet opening profile /, a milder 

deceleration g is needed to produce an increase in a IX than to 
separate the X-ray and radio peaks. This is in agreement with 
observations that show that the decline of the radio-to-X-ray 
flux is a more widely observed trend that the X-ray-radio peak 
displacement. 

For the rest of this work, we consider conical jets with 
/ = (1 + az)/( 1 + az 0 ), where a = 0 corresponds to a cylin- 
drical flow, and deceleration profiles g = (z/z 0 )~ e , where e > q s 
0. The behavior of the EC and synchrotron flux profiles also 
depends on the observing angle. To demonstrate this, we plot 
in the upper (lower) panel of Figure 1 the EC (synchrotron) 
flux along a decelerating jet, for a range of angles. It can be 
seen that while for small angles the EC flux peaks at the base 
of the jet, at larger angles the peak shifts downstream as a result 
of the fact that the emission from the faster base of the jet is 
beamed out of the line of sight [for a given 6 , A(z) peaks at 
T(z) ~ 1/ sin 6}. The behavior of the synchrotron flux is also 
interesting. At angles 6 ^ 0 the peak of the synchrotron emis- 
sion shifts downstream, and for a given angle it is farther down- 
stream than the EC peak, in agreement with observations. 

It is interesting to see how much the above considerations 
change when radiative losses are important. The simplest case 
is that of a perfectly collimated (i.e., ot = 0) cylindrical jet. 

This is also the case that requires the mildest deceleration for 
producing a given a rx increase or X-ray-radio separation. As- 
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Fig. 2. — Radio (top), optical (middle). X-ray (bottom) maps and the radio- 
to-X-ray spectral index (below the X-ray map) for a decelerating flow 
observed under 0 = 6 ° (left) and 0=12° (right). The cylindrical flow, assumed 
to be at Z = 0.3, decelerates from T 0 = 6 to T = 2 with e = 1 and B 0 = 
3 x 10 -6 as it propagates a distance of 200 kpc. The EED is a power law of 
index 5 = 2 and 7 ^ = 10 6 at the base of the flow. This calculation includes 
radiative losses. Each map is normalized, with red indicating the maximum 
luminosity per unit area. The projected jet lengths are *21 kpc for 0=6° and 
=s42 kpc for 0 = 12 °. 


suming an EED slope 5 = 2, the EED remains a power law 
of the same slope following equation (1), with a maximum 
energy electron Lorentz factor 


7max(Z) 1 + Y ra „ (I JC 1 (z" 0,+3V3 - 1) + C 2 (z * -2 ‘ v3 - D] ’ 


(4) 


where y maxf , 0> is the upper cutoff of the EED at z 0 , z* = z/z 0? 
C 1 = 3z o ABj/8ir(10e + 3), C 2 = 4z 0 AU(l + Z) 4 T 0 2 /( 3 - 2e), 
and A = 4a T /3m e c 2 T 0 . The maximum EC energy e EC max as a 
function of z is then e EC max = e 0 A 2 (z)y 2 zx (z), while the max- 
imum synchrotron energy is e 5 , max = B 0 z*A(z)yl„(z)/B cr For 
a given observing energy and observing angle, one can find 
the radiating part of the jet and perform the optically thin 
radiative transfer by integrating equations (2) and (3) along 
each line of sight for a two-dimensional array of lines of sight 
that covers the source as projected on the observer’s plane of 
the sky. 

The results of such calculations are shown in Figure 2, w here 
it can be seen that the optical emission is confined at the base 
of the flow because of the strong radiative losses of the high- 
energy electrons, practically marking the site of strong particle 
q 6 acceleration. The morphology of the radio emission is angle- 
dependent: at 0 = 6° it covers the entire extent of the jet, while 
at 0 = 12° the emission of its inner part is dimmed because 
most of the radiation is beamed outside our line of sight. The 
X-ray emission at both angles peaks close to the base of the flow 
because of the decrease of the comoving CMB photon energy 
density away from the flow base. Regarding the relative radio- 
X-ray morphology, while at 0 = 12° the offset betw een the tw'o 
images is very large, at 6 = 6° the tw'o images overlap, although 
the radio image clearly extends farther downstream. Also, as the 
radio-X-ray spectral index plots show, in both cases ol tx increases 
downstream in agreement with observations. 

We now compare this picture with the one we w'ould get 
from the same system if radiative losses were balanced by 
distributed reacceleration (Fig. 3). The X-ray images remain 
practically the same because the radiative losses for the EC X- 
ray-emitting electrons in this Z = 0.3, T 0 = 6 jet are negli- 
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Fig. 3. — Same as in Fig. 2, but with radiative losses balanced by distributed 
reacceleration and without the optical map, which in this case would be iden- 
tical to the radio map. 
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gible. The radio image, however, changes significantly, partic- 
ularly so for the 6 = 6° orientation. Now that the radio-emitting 
electrons retain their energy, the outer part of the jet is much 
brighter than the inner part, resulting in a clear separation be- 
tween the radio and X-ray images. 

3. DISCUSSION AND CONCLUSIONS 

We have proposed that the increase of the radio-to-X-ray 
flux ratio along the length of the jets of powerful quasars and 
the occasional offset of the jet images in these wavelengths are 
naturally accounted for in terms of relativistic flows that de- 
celerate over the entire length of the jet. Despite this decel- 
eration, the jets remain relativistic (T of about a few) to their 
terminal hot spots (Georganopoulos & Kazanas 2003), within 
which they eventually attain subrelativistic speeds. Our pro- 
posal provides, for the first time, a means for deducing the jet 
kinematics through simple models of their multiwavelength 
images; these can then be checked for consistency when cou- 
pled w'ith the detailed hot spot emission, which, as shown in 
earlier w'ork (Georganopoulos & Kazanas 2003), depends on 
the value of T at this location. We note here that a change of 
the Doppler factor A can also result from a jet curving aw'ay 
from (or into) the observer’s line of sight without actually 
decelerating. In this case, due to the stronger dependence of 
ECjon A (Dermer 1995; Georganopoulos et al. 2001), it is also q 7 
possible to obtain increasing (or decreasing, which is not ob- 
served) radio-to-X-ray ratios; however, because the comoving 
CMB photon density and magnetic fields remain constant in 
this case, such an arrangement cannot produce the observed 
displacements between the X-ray and radio peak emissions. 

The scenario we propose here finds further support from 
modeling jets with several sequential individual knots: Sam- 
bruna et al. (2001) noted the need for a gradual decrease of 
the Doppler factor and/or an increase of the magnetic field in 
order to reproduce the emission from the knots along the jet 
of 3C 273 with simple one-zone models. These knot-to-knot 
variations can be naturally incorporated within the context of 
a decelerating collimated flow, as w'e propose. y 

These same maps also may or may not indicate the need for ^ — 
distributed particle reacceleration along these jets, when the 
EC loss length scale of their radio-emitting electrons [oc(l 4- q s 
Z)- 4 r~ 2 ] j s shorter than the observed length of the radio jet. 

This appears to be the case with sources like PKS 1127 — 145 
at Z = 1.187 (Siemiginow'ska et al. 2002) and possibly GB 
1508+5714 at Z = 4.3 (Yuan et al. 2003; Siemiginowska et 
al. 2003; Cheung 2004), which show' their peak radio emission 
displaced from that of the X-rays: in the absence of reacce- 
leration, EC losses of the radio-emitting electrons w'ould pro- 






PROOF 4 


GEORGANOPOULOS & KAZANAS 


Vol. 604 


duce radio jets shorter than X-ray jets. Much about the reac- 
celeration process of these jets can be inferred from their 
morphology at different frequencies. For example, the knotty 
optical jet morphology shows that reacceleration is not strong 
enough to offset the EC-dominated losses of the optically emit- 
ting synchrotron electrons (y ~ 10 6 -10 7 ). 

The radiative efficiency of these jets is less than a few percent 
(e.g., Tavecchio et al. 2000), and the energy lost in deceleration 
must be either used to heat up the matter in the jet or transferred 
to material that is entrained by the jet. While the first option 
would result in an expansion of the jet, contrary to what is 
seen, entrainment from an external medium would load the jet 
with baryonic mass while decelerating it. Entrainment would 
produce velocity gradients across the jet, and in this sense a 
faster spine and a slower sheath are to be expected. However, 
the observed X-ray and radio jet morphologies suggest that the 
dominant effect must be a deceleration along the jet. A fast 
spine that does not decelerate substantially and carries most of 
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the jet power would produce a constant X-ray EC and radio 
synchrotron flux along the jet, in disagreement with observa- 
tions. A consequence of entrainment would be that even if the 
jet did not start as a baryonic one, entrainment would gradually 
enrich it with baryons and eventually a fraction AT/T of its 
power, where AT is the decrease of the Lorentz factor, would 
be carried by the entrained baryonic matter. This in turn would 
increase the radio lobe equipartition energy content estimates 
derived under the assumption of a leptonic composition. Our 
findings point to a picture in which powerful relativistic jets 
decelerate, depositing most of their power in their surroundings 
in the form of kinetic energy, their observed radiation being 
only the tip of the iceberg. 
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